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We present a theoretical description of Bernstein modes that arise as a result of the coupling be- 
tween plasmon-like collective excitations (upper-hybrid mode) and inter-Landau-level excitations, 
in graphene in a perpendicular magnetic field. These modes, which are apparent as avoided level 
crossings in the spectral function obtained in the random-phase approximation, are described to 
great accuracy in a phenomenological model. Bernstein modes, which may be measured in inelas- 
tic light-scattering experiments or in photo-conductivity spectroscopy, are a manifestation of the 
Coulomb interaction between the electrons and may be used for a high-precision measurement of 
the upper-hybrid mode at small non-zero wave vectors. 

PACS numbers: 78.30.Na, 73.43.Lp, 81.05.ue 



I. INTRODUCTION 

Most of the electronic properties of graphene, two- 
dimensional (2D) graphite, may be understood within 
the picture of non-interacting or weakly-correlated mass- 
less Dirac particles.^ A notable exception is the recently 
observed fractional quantum Hall effect,^ ? which arises 
in a strong magnetic field when a Landau level (LL) is 
only partially filled and when the Coulomb interaction 
becomes the relevant energy scale due to a quenched ki- 
netic energy, similarly to the usual 2D electron gas in 
semiconductor heterostructures. Apart from this rather 
particular situation, the role of the Coulomb interactions 
may be quantified with the help of the graphene fine- 
structure constant aa = e^/hevp ~ 2.2/e, in terms of 
the Fermi velocity vp ~ 10® m/s and the dielectric con- 
stant e of the surrounding medium. 

Even if the value of ac indicates that graphene should 
be a moderately correlated material, experimental indica- 
tions for the role of Coulomb interactions in the absence 
of a magnetic field are rather sparse. Recent angular- 
resolved photoemission experiments have revealed fea- 
tures that hint at "plasmaron" excitations that may be 
viewed as quasi-particles dressed by interaction-induced 
plasmon excitations.^ Furthermore, an indirect determi- 
nation of the screened fine structure constant, defined 
as aQ(q, w) — e'^ /hvpe{q,uj), where the dynamic screen- 
ing due to inter-band processes is encoded into the di- 
electric function e(q, w), has been measured by means of 
X-ray scattering in graphite. In the w — >■ and long 
wavelength limit, the value given in Ref. "Ffor graphene 
(extracting indirectly the polarization of graphene from 
measurements on graphite) is ajj ~ 0.14, i.e. a value that 
is roughly one order of magnitude smaller than the one 
mentioned above, for the case of freestanding graphene 
(e ~ 1). The investigation of interaction-related effects 
may therefore be viewed as one of the major issues in 
fundamental research on graphene. 

Interaction effects are expected to play a role in 
graphene exposed to a strong perpendicular magnetic 



field in the form of plasmon-like excitations.'^'l^ Most 
saliently, the relativistic character of the electrons in 
graphene gives rise to rather exotic linear magncto- 
plasmons,'^in addition to the upper-hybrid mode (UHM) 
which is a mixture of the cyclotron mode and the usual 
2D plasmon with its low-energy ^ dispersion, in terms 
of the wave vector q. In the context of the 2D elec- 
tron gas in semiconductor heterostructures, the UHM is 
coupled via the Coulomb interaction to inter-LL excita- 
tions, and their hybridization gives rise to avoided level 
crossings, known as Bernstein modes.'i^ Bernstein modes 
have been measured in photo-conductivi ty spe ctroscopy 
and inelastic light-scattering experiments, 2^"^^ and these 
techniques may in principle also be applied to graphene. 

In this paper, we study theoretically the Bernstein 
modes in graphene. These Bernstein modes become ap- 
parent in the spectral function, which takes into account 
the Coulomb interaction between the electrons in the 
random-phase approximation (RPA) . We propose a phe- 
nomenological model that captures the relevant features, 
as the position and size of the anticrossings between the 
UHM and the inter-LL excitations, that appear in the ex- 
citation spectrum of this system. Moreover, our results 
could be used to analyze future experimental results and 
to identify the different modes. For this aim, we propose 
two different experimental setups that could be applied 
to graphene in order to measure the Bernstein modes as 
well as the dispersion relation of the UHM in this mate- 
rial. In contrast to 2D electrons in semiconductor het- 
erostructures with a parabolic band dispersion and an 
equidistant LL spacing, one observes a plethora of inter- 
LL transitions in graphene as a consequence of the \^/ Bn 
scaling of the relativistic LLs, where A = + or — for the 
conduction and the valence band, respectively, and the 
integer n denotes the LL index. This results in a large 
number of intersections between the inter-LL transitions 
and the UHM, such that a measurement of the Bernstein 
modes would in principle allow for a high-precision deter- 
mination of the dispersion relation of the upper-hybrid 
mode. 
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FIG. 1. (Color online) RPA Spectral function for electrons in 
graphene, with Nf = 3, 6 = Q.Q2vf /Ib, for a fine-structure 
constant ac ~ 1- The dashed black lines are obtained from 
the expression ([8| for the coupled modes, for the same param- 
eters and 7 — 3/4. The continuous lines delimit the region 
for intra-band (I) and that for inter-band (II) excitations. 



The paper is organized as follows. In Sec. [TTj we dis- 
cuss the Bernstein modes which are visible in the RPA 



spectral function (Sec. II A) and propose a phenomeno- 
logical model (Sec. II B) that allows for a quantitative 



description of the position and the strength of the modes. 
Section |III| is devoted to possible experimental observa- 
tions of Bernstein modes in graphene, and we present our 
conclusions in Sec. IIVI 



phenomenological level. The form factors J^\n,\'n' (q) are 
due to the graphene-LL wave functions and take into ac- 
count the chirality of the carriers. They may be expressed 
in terms of the cyclotron variable r] — r — R, where r is 
the position of the electron and R that of the center of 
the cyclotron motion, which is a constant of motion, 

^An,A'«'(q) = i;!,;- {n - 1 le-^i-^^l n' - l) 

+XXXK'{n\e''''-'^\n'), (3) 



and 

v/(l - 5n,o)/2 and 2: 



we have used the short-hand notation 1* = 
= ^{l + Snfi)/2. The states |n) 



are the usual eigenstates of the harmonic-oscillator oper- 
ators, a'a\n) = n\n), where a = {rj^ + iriy)/^/2lB- 

The spectral function, i.e. the imaginary part of the 
RPA polarizability n"(aj, q) /en.pA{^, q), is shown in Fig. 
[T] where we have chosen, for illustration reasons, a value 
of ac = 1 for the graphene fine-structure constant. The 
Fermi level Ep = Xw'y/Wp is fixed in the LL Np = 3 
in the conduction band (A = +), and we have chosen an 
impurity broadening oi S — 0.02vp/Ib- The main spec- 
tral weight is concentrated in the UHM, which disperses 
aJl3 



(4) 



and which may be viewed as the high-field descendent of 
the B = 2D plasmon with the approximate dispersion 
relatior 



II. THEORETICAL DESCRIPTION OF 
BERNSTEIN MODES IN GRAPHENE 

A. Bernstein modes in the RPA spectral function 

Quite generally, the collective charge excitations gen- 
erated by the Coulomb interaction between the electrons 
may be obtained from the zeros of the RPA dielectric 
function 
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eflPA(w,q) = 1 - -^n°(a;,q). 



(1) 



in terms of the polarizability 11° (w, q) for non-interacting 
2D electrons. For electrons in graphene in the integer 
quantum Hall regime, the pola rizability reads (we use a 
system of units with h = l|2El 



n°(u;,q)= 



X 



g[np{Xn) — np{X'n'y\ 
u!'{Xy/n - X'Vnf ) +Ljj + i5 

(q)P 



(2) 



where np{Xn) is the Fermi distribution function, w' = 
^/2vp/Ib is the characteristic LL frequency, g — ^ rep- 
resents the four-fold spin-valley degeneracy, and 5 takes 
into account the disorder-induced LL broadening on a 



. , 2e'^Ep 

^piq) ^ \l — — q 



iv^^q^. 



(5) 



The first term in this expression yields the usual ^Jq dis- 
persion of the classical plasmon,!^ which is clearly visible 
in Fig. [T] while the second one represents a higher-order 
correction due to quantum effects. The small-g expan- 
sion of the RPA polarization function yields a negative 
prefactor for the quantum corrections, 7 — — a^l^Sl^ in 
contrast to non-relativistic electrons. It seems, however, 
that 7 depends itself on q and ac namely at larger values 
of the wave vector where 7 crosses over to positive values. 
Here, we use 7 as a fitting parameter to our RPA results 
in Fig. [1] where the UHM is best described by the value 
7 = 3/4 as noted by Shung,!^ in the shown wave- vector 
range and for an interaction parameter of ac = 1. The 
main effect of the magnetic field is the q = Q gap, which 
the plasmon acquires and that is given by the density- 
dependent cyclotron gap ujc{B) = eBvp/Ep in Eq. m. 

It is clearly visible in Fig. [TJthat, although the main 
part of the spectral weight is found in the upper-hybrid 
mode, this mode is not a continuous line. Instead, one 
notices a series of avoided level crossings whenever the 
upper-hybrid mode coincides with the energy of a dis- 
persionless inter-LL transition, riAn.n' — lo' {\fn! — X\fn)^ 
where an electron is promoted from the LL n in the band 
A to n' in the conduction band (above the Fermi level). 
These avoided level crossings are nothing other than the 
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Bernstein modes, which occur at the wave vectors 

qls = -V2N^ 

1 



(6) 



"0 + 7 



Nf 



47V| 



in graphene. 



B. Phenomenological model 

In order to describe the couphng between the UHM 
and the inter-LL excitations within a phenomenological 
model, we treat the UHM as a bosonic excitation, de- 
scribed by the Hamiltonian 



Huh = y^^ufe(g, g)6|,6q. 



Here, the boson operator &q' is proportional to the 
Fourier component Puh{^) of the density operator that 
constitutes the plasmon-type mode. This mode is cou- 
pled via the Coulomb interaction 



H. 



coupi = 7 n~T + Pu'i(-q)p(q)] 

4 ^-^ e|q| 



to the inter-LL excitations described by the density com- 
ponents 



(q) V(m|e-^'i-^|m')cl„_^cv„'.™', 



Xn.X'n' 



in terms of the fermionic electron operators c^^2 m- This 
means that the electronic density is separated into a part 
that forms the UHM and another one that describes the 
inter-LL excitations, and formally one needs to introduce 
a constraint to avoid double counting of the electronic 
degrees of freedom. 

The Coulomb coupling between the UHM and the LL 
excitations yields a renormalization of ujuh that may be 
calculated via a Dyson equation for the dressed propa- 
gator of the UHM, similarly to the magneto-phonon res- 
onance discussed in Ref. ^22, The avoided level crossing 
is then governed by the polarizability ([2]), which is dom- 
inated by the resonant term with uj ~ ujuh — ^\n,n'^ and 
the equation giving the poles of the dressed propagator 
reduces to 
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(7) 



where V = {e'^ /eql%)\J^\n.7i'{(l)\'^ is the effective coupling 
constant. From Eq. ([?]), one obtains the two solutions 



(8) 
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The solutions are plotted in Fig. [T] in the form of 
dashed black lines, and one notices the good agreement 
with the maxima of spectral weight obtained from Eq. 

At resonance [tUuhiq) — f^An,ri'], Eq. ^ may be 
linearized. 



UJ± 



in terms of the relative splitting parameter 



'(<?) = 



^/gV 

vf/Ib 



qls 



(q)l' 



(9) 



One notices in Fig. [T] that the UHM couples more 
strongly to inter-LL excitations that emerge from the 
inter-band regime (region H), with A = — , than to those 
associated with intra-band transitions (region I). Indeed, 
the relative splitting parameter ([9| is proportional to the 
|-^An,n'(q)P: which is important in the regions of the 
spectrum which correspond to the particle -hol e contin- 
uum in the absence of the magnetic field. I^lifll Because 
the UHM does not enter the intra-band region (I) of this 
continuum, the coupling to the corresponding inter-LL 
excitations (with A = +) is strongly suppressed, in con- 
trast to those in the inter-band part of the particle-hole 
continuum, where the upper-hybrid mode enters. 



III. DISCUSSION AND IMPLICATIONS FOR 
EXPERIMENTS 

Experimentally, Bernstein modes have been ob- 
served in 2D electron systems in semiconductor 
heterostructures.'^^Hm in contrast to three-dimensional 
metals,'22l the main difficulty stems from the vanishing 
plasmon dispersion at q = 0, such that a direct measure- 
ment of Bernstein modes in homogeneous 2D systems via 
spectroscopic means is impossible. It is therefore neces- 
sary to impose a non-zero value of the wave vector to 
the system. Two different techniques have been success- 
fully used to do so. The first one consists of a grated 
coupler, with a well-defined periodicitya, parallel and in 
close vicinity to the 2D electron gas.'i2"i^The wave vector 
go = 27r/a in the 2D electron gas, at which one may mea- 
sure the electromagnetic response, is thus fixed by the 
periodicity of the grid via an electromagnetic coupling 
effect. This allows one to measure collective charge exci- 
tations via transmission spectroscopy (in the far-infrared 
regime) at a characteristic (non-zero) value of q. 

An alternative technique consists of inelastic light- 
scattering experiments at a well-defined angle 9 between 
the incident photon and the vector perpendicular to the 
2D electron gas.^ In a light-scattering process, the pho- 
ton thus transfers a momentum gg — Awsin^/c to the 
electrons, where Aw = {cui — uJs) is the energy difference 
between the incident and the scattered photon, and c is 
the velocity of light. A resonance at the wave vector q 
and energy Aw is then interpreted as a collective charge 
excitation, such as a plasmon. 
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FIG. 2. Bernstein modes as a function of the magnetic field 
obtained from Eq. ([8|, for a fixed value oi qo/kp = 1 and 
aa = 1- Dashed lines represent inter-LL excitations, and the 
continous line shows the UHM. The inset shows a zoom on 
the high-field region of the avoided level crossings. 



In both techniques, which may in principle also be ap- 
plied to graphene, it turns out to be simpler to study the 
Bernstein modes at the resonance condition ([6]) as a func- 
tion of the magnetic field while keeping the wave vector 
fixed. One notices two essential differences of graphene 
with respect to a non-relativistic 2D electron gas in semi- 
conductor heterostructures. First, because the electrons 
in graphene reside at the surface, the system has the 
advantage of being directly accessible by surface spec- 
troscopic means, as e.g. scanning tunneling microscopy. 
Second, as a consequence of the non-equidistant LL spac- 
ing in graphene, the inter-LL transitions are not multiples 
of the fundamental cyclotron frequency uc = eB /mB, in 
terms of the band mass mB of the semiconducting host 
material. Indeed, in a 2DEG there are, in general, many 
electron-hole transitions contributing to the same energy, 
due to the equidistant LLs. However, each particle-hole 
process contributes to a different energy in graphene be- 
cause of the V Bn dispersion of graphene LLs. As a conse- 
quence, within the same energy window, there is a higher 
number of anticrossings between the UHM and the inter- 
LL transitions (Bernstein modes) in graphene as com- 
pared to a 2DEG.^'^ Because the UHM is determined by 
the avoided level crossings with intcr-LL excitations, this 
would, in principle, allow for a more precise measurement 
of the UHM in graphene as compared to non-relativistic 
electron systems. 

In Fig. [2j we have plotted the solutions ([s]) of the Bern- 
stein modes as a function of the magnetic field instead 
of the wave vector in order to make a closer connection 
with possible experiments. The wave vector is fixed at 
qo/kp = 1 and we have chosen again aa = 1 for illus- 
tration reasons. Notice that, in contrast to the above 
discussion, the natural units are no longer vf/Ib for 
the energy and 1^^ for the wave vector but the density- 



dependent Fermi energy Ep 
vector kp — y/Awn^afg 



vpkp and Fermi wave 
y/Tm^i. Both units are there- 
fore fixed by the carrier density riei, whereas the in- 



dex Np changes as a function of the magnetic field as 
Np ~ k\'l%/2, where we have neglected the discrete- 
ness of Np. The inter-LL excitations now disperse as 
\/B{\/n-\-\/n')^ whereas the UHM is only weakly dispers- 
ing in the plotted magnetic-field range. As one expects 
from Eq. (|9|, the Bernstein modes are most prominent 
at high magnetic fields, and the position of the avoided 
level crossings are easily obtained from Eq. ([6]). Neglect- 
ing the quantum corrections to the plasmon frequency 
(7 — 0) in the small wave-vector limit {qo/kp ^ '^ctc/l), 
the Bernstein modes are expected at the field values 



(eB)- 



{Vn' - x^y 

2aGqokp 



(10) 




2^/2. 



acqokp 



Notice that the zero-field limit of the dispersion of 
the UHM, which one may extrapolate from the mea- 
surement of the Bernstein modes, gives direct access 
to the bare graphene fine-structure constant aa- In- 
deed, if one expresses Eq. ^ in terms of fc^ and ac, 
[LOuh{q,B — > 0)/vpy — 2aQkpq, one may determine the 
fine-structure constant in the small-g limit from the ex- 
pression 



1 



ul^{qo,B^O) 



2y^'Kneiqo 



(11) 



However, we emphasize that it is the unscreened 
graphene fine-structure constant which occurs in these 
expressions and not the screened one a^. Naturally, the 
information about aa is also encoded in the size of the 
mode splitting ([9]), but it seems more delicate to extract 
because of the complicated B-field and q dependence of 
the form factor. 



IV. CONCLUSIONS 

In conclusion, we have theoretically analyzed the Bern- 
stein modes in graphene, which arise as a consequence of 
the interaction between plasmon-type modes and inter- 
LL excitations. The proposed model captures the posi- 
tion and strength of the couplings, as compared to the 
numerical solution of the RPA polarizability. Finally, we 
have shown that a possible measurement of Bernstein 
modes with the help of techniques that have been suc- 
cessfully used in conventional 2D electron systems with a 
parabolic band dispersion^^^^^ and that may equally be 
used in graphene would allow for a measurement of the 
dispersion relation of the UHM. 
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